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One method to reconstruct the scalar field potential of infiation is a 
perturbative approach, where the values of the potential and its derivatives 
are calculated as an expansion in departures from the slow-roll approxima- 
tion. They can then be expressed in terms of observable quantities, such 
as the square of the ratio of the gravitational wave amplitude to the den- 
sity perturbation amplitude, the deviation of the spectral index from the 
Harrison-Zel'dovich value, etc. Here, we calculate complete expressions 
for the second-order contributions to the coefficients of the expansion by 
including for the first time corrections to the standard expressions for the 
perturbation spectra. As well as offering an improved result, these cor- 
rections indicate the expected accuracy of the reconstruction. Typically 
the corrections are only a few percent. 
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I. INTRODUCTION 



An intriguing prospect raised by recent large-scale structure observations, par- 
ticularly those of the Cosmic Background Explorer (COBE) satellite is that ob- 
servations may soon provide rather detailed information regarding the nature of the 
vacuum energy driving inflation 0,0,^. In most models this vacuum energy is iden- 
tified as the self-interaction potential of a scalar inflaton field, and its precise form is 
determined by some particle physics model. Since there currently exist many possible 
models 0, it is of great interest to investigate whether observations can select which, 
if any, of these models is correct. As well as having detailed implications for the 
initial conditions for structure formation in the Universe, the energy scale of inflation 
provides a link with particle physics at high energies, and may be a useful method of 
probing models of unification. 

The prime observational consequences of inflation derive from the stochastic spec- 
tra of density (scalar) perturbations and gravitational wave (tensor) modes generated 
during inflation. Each stretches from scales of order centimeters to scales well in ex- 
cess of the size of the presently observable Universe. Once within the Hubble radius, 
gravitational waves redshift away and so their main influence is on the large-scale 
microwave background anisotropies, such as those probed by COBE 0. Advanced 
gravitational wave detectors such as the proposed beam-in-space experiments may be 
able to detect the gravitational waves on a much shorter (about 10^'^cm) wavelength 
range |]^. The density perturbations are thought to lead to structure formation in 
the Universe. They produce microwave background anisotropies across a much wider 
range of angular scales than do the tensor modes, and constraints on the scalar spec- 
trum are also available from the clustering of galaxies and galaxy clusters, peculiar 
velocity flows and a host of other measurable quantities . 

Recently, we provided a formalism which allows one to reconstruct the inflaton 
potential ^(0) directly from a knowledge of these spectra This developed an 

original but incomplete analysis by Hodges and Blumenthal |T0|. An important result 
that follows from our formalism is that knowledge of the scalar spectrum alone is in- 
sufficient for a unique reconstruction. Reconstruction from only the scalar spectrum 
leaves an arbitrary integration constant, and since the reconstruction is nonlinear, dif- 
ferent choices of this constant lead to different functional forms for the potential. A 
minimal knowledge of the tensor spectrum, say its amplitude at a single wavelength, 
is sufficient to lift this degeneracy. With further information the problem becomes 
overdetermined, providing powerful consistency relations which would exclude infla- 
tion if not satisfied. 

The most ambitious aim of reconstruction is to employ observational data to 
deduce the complete functional form of the inflaton potential over the range cor- 
responding to large-scale structure. The observational situation is some way from 
providing the quality of data that this would require, and at present a more realistic 
approach is to attempt a reconstruction of the potential about a single point 0o H. 
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For this one requires such information as the amphtudes of both scalar and tensor 
modes and also the spectral index of the scalar perturbations at a single scale. It 
is possible that such information can be deduced from a combination of microwave 



background experiments that span a range of angular scales [|lT|. In this paper it is 



our aim to provide an improved calculation of the coefficients of such a perturbative 
reconstruction. 

To some extent all inflationary calculations rely on the use of the slow-roll approx- 
imation. In the form we present here, the slow-roll approximation is an expansion in 
terms of quantities deflned from derivatives of the Hubble parameter H. In general 
there are an inflnite hierarchy of these which can in principle all enter at the same 
order in an expansion. However, to calculate V{(j)Q) one needs only the flrst deriva- 
tive of if, for V\(Pq) one needs up to the second derivative of H, and so on. These 
parameters can be converted into more observationally related quantities, as we shall 
see. 

The slow-roll approximation arises in two separate places. The flrst is in simpli- 
fying the classical inflationary dynamics of expansion, and the lowest-order approxi- 
mation ignores the contribution of the inflaton's kinetic energy to the expansion rate. 
The second is in the calculation of the perturbation spectra; the standard expres- 
sions are true only to lowest-order in slow-roll. In our earlier work P,^, we utilized 
the Hamilton- Jacobi approach |12| to treat the dynamical evolution exactly, but were 
forced for analytic tractability to retain the lowest-order approximation for the pertur- 
bation spectra. Technically therefore, the results were accurate only to lowest-order, 
though in models close to the power-law inflation limit this hybrid approach offers 
substantial improvements for certain quantities. 

Until recently further improvements have not been possible, but a very elegant cal- 
culation of the perturbation spectra to next order in slow-roll has now been provided 
by Stewart and Lyth [|I^ . This does not permit analytic progress in functional recon- 
struction, but their results can be combined with the Hamilton- Jacobi approach to 
generate the complete second-order term in perturbative reconstruction. The purpose 
of this work is to calculate this correction. This serves two useful purposes. Firstly, 
the results allow a more accurate reconstruction to be performed, and secondly the 
relative size of lowest-order and second-order contributions provides a useful (though 
not rigorous) measure of the theoretical error in reconstruction. As we shall see, even 
the lowest-order results are typically accurate to within a few percent. 



II. TO SECOND-ORDER IN SLOW-ROLL 

We shall use the notation and philosophy of our earlier paper 0, except that the 
perturbation spectra shall be given by expressions improving on Eq. (3.4) of that 
work. The Hamilton- Jacobi equations arise when one uses the scalar fleld as a time 
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variable, and writes the Hubble parameter H = a/a, where a is the scale factor, as a 
function of 



The field equations are |I2 
3 
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-\K'V{ct>) 

-2H' , 



(2.1) 
(2.2) 



where dots are time derivatives, primes are derivatives, = Svr/mp; and mpi is the 
Planck mass. Without loss of generality we may assume > 0, so that H'{(j)) < 0. 
Where square roots appear later this choice is used to fix the sign of the prefactor. 

The slow-roll approximation can be specified by parameters defined from deriva- 
tives of H{(j)). There are in general an infinite number of these as each derivative 
is independent, but usually only the first few enter into any expressions. We shall 
require the first three, which are all of the same order when defined byQ 
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The slow-roll approximation applies when these slow-roll parameters are small in 
comparison to unity. The condition for infiation, a > 0, is precisely equivalent to 
e < 1. 

The lowest-order expressions for the scalar {As) and tensor {Aq) amplitudes as- 
sume {e, ^} are negligible compared to unity. Improved expressions for the scalar 
and tensor amplitudes for finite but small {e, r/, ^} were found by Stewart and Lyth 



A. 



A 



G 



//[l-(C + l)e] , 



47^3/2 



(2.4) 
(2.5) 



where C = —2 + ln2 + 7 ~ —0.73 is a numerical constant, 7 ^ 0.577 being the Euler 
constant. The right hand sides of these expressions are evaluated when the scale in 
question crosses the Hubble radius during infiation, 27r/A = aH. The spectra can 
equally well be considered to be functions of wavelength or of the scalar field value. 
Eq. ( |2.7| ) below allows one to move from one to the other. 

The standard results to lowest-order are given by setting the square brackets to 
unity. Historically it has been common even for this result to be written as only 



^Let us stress that our choice > imphes ^/e = —^/2/f?H'/H: one needs to be careful with 
the signs to reproduce our results. 
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an approximate equality (the ambiguity arising primarily because of a vagueness in 
defining the precise meaning of the density perturbation), though the precise nor- 
malization to lowest-order was established some time ago by Lyth (see also the 
discussion in [Q). 

The improved expressions for the spectra in Eqs. ( p.4| ) and ( |2.5| ) are accurate in 
so far as e and r] are sufficiently slowly varying functions that they can be treated 
adiabatically as constants while a given scale crosses outside the Hubble radius. Cor- 
rections to this would enter at next order. This differs from the usual situation in 
which H is treated adiabatically. For the standard calculation to be strictly valid H 
must be constant, but provided it varies sufficiently slowly (characterized by small e 
and I77I), it can be evaluated separately at each epoch. This injects a scale depen- 
dence into the spectra. There is a special case corresponding to power-law inflation 
for which e and rj are precisely constant and equal to each other. In this case the 



above expressions for the perturbation spectra are exact |T^JT^. Furthermore, the 
corrections to each spectrum are the same and they cancel when the ratio is taken. 
In the general case e and t] may be treated as different constants if it is assumed that 
the timescale for their evolution is much longer than the timescale for perturbations 
to be imprinted on a given scale. This assumption worsens as rj is removed from e, 
which would be characterized by the next order terms becoming large. 

Throughout we shall be quoting results which feature a leading term and a cor- 
rection term linear in the slow-roll parameters. We shall utilize the symbol "~" 
to indicate this level of accuracy throughout. The correction terms shall be placed 
in square brackets, so the lowest-order equations can always be obtained by setting 
the square brackets equal to one. A useful relationship can be obtained from Eqs. 



e ^ 41 [1 - 2C(e - r^)] . (2.6) 

As we shall see, rj is encoded in the spectral index of the scalar perturbations, whose 
deviation from unity must also be small for slow-roll to apply. 

A key equation in Refs. is the consistency equation, which connects the scalar 
spectrum to the tensor spectrum and its derivative. The spectra as given in Eqs. (|2.4|) 



and (|2.5| ) are functions of the value of when the fluctuations crossed the Hubble 
radius during inflation. This is converted into a dependence on wavelength A with 
the relation ^ 

Differentiation of Eq. (|2.5|) with respect to implies that 
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and it follows that 
A dAc A% 



[1 + 3e - 2r]] . (2.9) 



Ac d\ A| 

As expected this agrees with the expansion of the corresponding expression in Ref. 
0] for the special case e = rj. This equation is interesting in its own right, but for 
perturbative reconstruction its use is restricted to the removal of derivatives of the 
tensor spectrum. 



III. PERTURBATIVE RECONSTRUCTION TO SECOND ORDER 



The aim now is to obtain expressions for the potential and its derivatives about 
a single point (po, given information regarding the spectra at the scale Aq which left 
the horizon at = 0o- The four main quantities of observational interest are the am- 
plitudes and spectral indices of the two spectra. However, in view of the consistency 
equation, Eq. (|2.9| ), only three of these are independent. We shall concentrate on the 
two amplitudes and the scalar spectral index, since these are probably the easiest to 
measure. The scalar spectral index n is defined by 



n 



d\nAj.{X) 
dlnX 



To lowest-order in slow-roll one can show that the spectral index is given by 



n 



(3.1) 



(3.2) 



which provides the route to determining rj. Conceptually we are passing from these 
three observables to the three parameters that describe the potential, which are the 
slow-roll parameters e and r/, and the overall normalization. In terms of the observ- 
ables, therefore, the slow-roll approximation amounts to an expansion in both A^/A^g 
and (1 — n), which give corrections to the same order. We shall see that the correction 
term for l^"(0o) requires the introduction of the third slow-roll parameter ^, requiring 
a new independent observable to determine it. 

One obtains directly from the field equation Eq. ( |2.1| ) and the definitions of the 
spectra in Eqs. (p.4|) and (p.5|) an expression for the amplitude of the potential: 



^(0O) 



487r' 
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l+[- + 2C 



A|(Ao) 



1 + 0.21 



(3.3) 
(3.4) 



In Refs. we gave the numerical factor on the second-order term as —1/3, which 
incorporated only the dynamical slow- roll corrections. In fact, the spectral corrections 
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to V{(j)o) dominate the dynamical ones for any inflation model, reversing the sign 
of the correction, which may be significant if the tensors are important. However 
the relative contribution of the scalar and tensor modes to large angle microwave 
anisotropics with our spectral normalization is given approximately by -R 0, where 



R 



2A| 
25AI 



(3.5) 



Even if the contributions to the anisotropics from scalar and tensor modes are equal, 
the correction term in the potential is only 2%. This is a powerful indication that 
even the lowest-order perturbative reconstruction promises to be very accurate. 

To obtain V'^cpo) we need the scalar spectral index at Aq, denoted by uq. Differ- 
entiation of the potential with respect to 0, followed by some straightforward algebra 
gives 



V'i 
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[l + (C + 2)e + (C-l/3)r/] , 
[1 



1.27e - IMr]] 
AUXo) , 



1-0.85 



0.53(1-%) 



(3.6) 



Note that for power-law inflation, which has 

25 AUXq) 
47r A^(Ao) 

the corrections in the square brackets nearly cancel, but other models |[T6|,p!7 



(3.7) 



can 



feature larger corrections, e.g. 16 % for natural inflation with no = 0.7. 

The calculation for y'(0o) is much more involved. One can show a precise rela- 
tionship 



V"( 



(3.8) 



A new observable will be needed to determine ^, the easiest example being the rate of 
change of the scalar spectral index. This would be substantially harder to measure, 
and it is fortunate that it only enters at second-order. [It would however enter at 
leading order in V"'{(j)o)]- From Eqs. ( ^.5] ) and ( |3.8| ), we can obtain the second-order 
correction to V"( 



A8tt- 



in terms of the slow-roll parameters 



■A'ciXo) (e + v) 



l + (2C + 2)e- 



3(e + r/) 



(3.9) 



It is however harder to convert the prefactor into observables, because there are now 
lowest-order terms in both e and f]. To generate the correct second-order term, it is 
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not enough to use the first-order expression for rj in terms of the spectral index. One 
must instead use the second-order result, as given by Stewart and Lyth |]13 



1 - n ^ 4e - 2r7 + 8(1 + C)e^ - (6 + 10C)er/ + 2Ce^ , 
(our 1] being the negative of their 6), which leads to 



(3.10) 



e + T] 



3e 



, 4C + 4 
1 + — z — e 



5C + 3 



1 + 
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4- 



-7] 



3^ 
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2C C-3 



3 3 
Substituting this into Eq. ( |3.9| ) yields 

1447r3 A^(A, 



1 - y^o 
2 ' 
1 - no 



247r3 



^ 4C + 10 C-3 

1 H e H f] + — f 

3 3 ' 3^. 



(3.11) 



(3.12) 



k2 



A^(Ao) (1 - no) [1 + (2C + 2)6] - ^A^^iv' + eO 



where the last term is entirely second-order. Note that there are two lowest-order 
terms. An interesting case is r] = — e, corresponding to oc 0^^^, for which the 
lowest-order term vanishes identically and the final term of Eq. ( |3.9| ) is the only one 
to contribute. The second derivative of the potential is the lowest derivative at which 
it is possible for the expected lowest-order term to vanish. 

The final step is to convert the second-order terms into the observables. As they 
are already second-order, one only needs the lowest term in their expansion to convert. 
From the expression for the spectral index, one finds to lowest-order that 



dn 



2e dlnX 



+ 5?7 - 4e . 



(3.13) 



Ao 



Note that the derivative of the spectral index is of order e^. To lowest-order we also 
have 



V 



2e 



;(1 - ^o); 



42 



(3.14) 



Progressively substituting all these into Eq. ( |3.12| ) yields 



--(1-no) 

3C - 1 dn 
^ 2 dlnX 



A|(Ao) 
- (12C-6) 



^(l-.o) + (36C + 2)||^ 



A - no) 



(3.15) 



Ao- 
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where the first two terms are lowest-order and the remainder are second-order.0 

For power-law models, the last term is zero and the remaining correction terms 
nearly cancel each other, though they are not small individually. For natural inflation 
models the correction terms are all individually small, with an overall correction of 
about 4% at hq = 0.7. (In natural inflation, (rfn/c? In A) |ao ~ r/^ ~ (1 — ^0)^/16 [p^ .) 



IV. DISCUSSION AND CONCLUSIONS 



To conclude, we have calculated the full second-order corrections to the perturba- 
tive reconstruction of the inflaton potential. The flrst-order terms agree with those 
we found previously P,^, while the second-order terms offer an improvement. They 
serve to quantify the expected errors in the perturbative reconstruction, and in general 
these errors are small. Even in cases where tensors provide a substantial contribu- 
tion to the large angle microwave background anisotropics and/or the spectral index 
deviates signiflcantly from unity, the corrections are typically only a few percent. 
Consequently, example flgures based on plausible data sets that we presented in our 
earlier papers remain valid. One then has some degree of confidence that one can use 



our lowest-order expressions as was done in |18 
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